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We present an analytic next-to-leading order QCD calculation of the partonic cross sections for 
single-inclusive lepton production in hadronic collisions, when the lepton originates from the decay 
of an intermediate electroweak boson and is produced at high transverse momentum. In partic¬ 
ular, we consider the case of incoming longitudinally polarized protons for which parity-violating 
single-spin asymmetries arise that are exploited in the W boson program at RHIC to constrain the 
proton’s helicity parton distributions. Our calculation enables a very fast and efficient numerical 
computation of the relevant spin asymmetries at RHIC, which is an important ingredient for the 
inclusion of RHIC data in a global analysis of nucleon helicity structure. We confirm the validity 
of our calculation by comparing with an existing code that treats the next-to-leading order cross 
sections entirely numerically by Monte-Carlo integration techniques. We also provide new compar¬ 
isons of the present RHIC data with results for some of the sets of polarized parton distributions 
available in the literature. 

PACS numbers: 12.38.Bx, 14.70.-e, 13.88.-|-e 


I. INTRODUCTION 


The W physics program at RHIC [1] is dedicated to 
providing new insights into the helicity structure of the 
proton. It exploits the violation of parity in the weak 
interactions, which gives rise to sm^Ze-longitudinal spin 
asymmetries in proton-proton collisions. The main fo¬ 
cus is on the production of W bosons, identified by 
their subsequent decay into a lepton pair. The charged 
lepton (or antilepton) is observed. From the corre¬ 
sponding cross sections for the various helicity settings 

(++), (H—), (—h), (-) of the two incoming protons one 

defines the spin asymmetry 



+ da^ — {da + da ) 
da++ + da+~ + da~+ + da 


dAa 

da 


( 1 ) 


As one can see, one takes the difference of cross sec¬ 
tions for positive and negative helicities of one proton, 
while summing over the polarizations of the other. The 
STAR collaboration at RHIC has published rather ex¬ 
tensive and precise data on last year [2], and new 

precise mid-rapidity data from PHENIX have just ap¬ 
peared [3]. Earlier measurements were reported by both 
PHENIX [4] and STAR [5]. Data sets with even higher 
statistics and kinematic coverage are expected in the near 
future. Typically, the data are presented at fixed rapidity 
of the charged lepton, which by convention is counted as 
positive in the forward direction of the polarized proton. 

It has long been recognized [6, 7] that offers ex¬ 
cellent sensitivity to the individual helicity parton distri¬ 
butions Alt, Am, Ad, Ad of the nucleon, where 

Af{x,Q^) = f+{x,Q^)-r{x,Q^), ( 2 ) 


with /■'■ (/“) denoting the distribution of parton / with 
positive (negative) helicity in a parent proton with pos¬ 
itive helicity. The distributions are functions of the lon¬ 
gitudinal momentum fraction x of the parton and of a 


“resolution” scale Q. Information on Am, Am, Ad, Ad 
is also accessible in (semi-inclusive) deep-inelastic lepton 
scattering (DIS) [8-12]. The key advantages of W boson 
production are that (i) it is characterized by momentum 
scales of the order of the W mass which are much higher 
than those presently relevant in DIS and hence deeper 
in the perturbative domain, (ii) it does not rely on the 
knowledge of hadronic fragmentation functions, thanks 
to its clean leptonic final state. In any case, information 
from the W program at RHIC is complementary to that 
from DIS. 

The main concept behind the RHIC measurements can 
be easily summarized: For W~ production, taking into 
consideration only the dominant ud -T W~ subprocess, 
the spin-dependent cross section in the numerator of the 
asymmetry in Eq. (1) is found to be proportional to the 
combination 

AM(a:i)d(a: 2 )(l — cosO)^ — Ad(xi)u(x 2 )(l + cos0)^ , (3) 

where for simplicity we have not written out the straight¬ 
forward convolutions over the parton momentum frac¬ 
tions. 0 is the polar angle of the negatively charged de¬ 
cay lepton in the partonic center-of-mass system, with 
d > 0 in the forward direction of the polarized parton. 
In the backward region of the lepton, one has X 2 
and 0 ^ 7r/2, so that the first term in Eq. (3) strongly 
dominates. Since the denominator of is proportional 
to u(xi)d(x 2 )(l — cos0)^ + d(xi)u(x 2 )(l + cos0)^, the 
asymmetry then provides a clean probe of Au(xi)/u(xi) 
at medium values of xi. By similar reasoning, in the for¬ 
ward lepton region the second term in Eq. (3) dominates, 
giving access to — Ad{xi)/d{xi) at relatively high xi. 

For production, within the same approximation, 
the spin-dependent cross section is proportional to 

Aj(a:i)M(x 2 )(l + cosd)^ — Am(xi)J(x 2 )(I — cosd)^ . (4) 

Here the distinction of the two contributions by con¬ 
sidering backward or forward lepton scattering angles is 




2 


less clear-cut than in the case of W~ because of the re¬ 
versal of the factors (1 ± cos0)^ relative to (3), which 
always suppresses the dominant combination of parton 
distributions. Therefore, both terms in (4) will com¬ 
pete. Nonetheless, the W~^ measurements at RHIC are 
of course of great value in the context of a global analysis 
of the helicity distributions. 

Given the importance of for constraining nu¬ 

cleon helicity structure, there has been a lot of activity 
on the calculation of higher-order QCD corrections to 
the relevant spin-dependent cross sections. Closed ana¬ 
lytic expressions for next-to-leading order (NLO) correc¬ 
tions to polarized W boson production were derived in 
Refs. [13, 14], with extensions to all-order resummations 
in [15, 16]. In these papers, direct observation of the W 
boson and its kinematics was assumed, which simplifies 
the calculation considerably but is not really applicable 
to the measurements at RHIC. The proper lepton de¬ 
cay kinematics was taken into account in three further 
studies [17-19]. The first two of these include the contri¬ 
butions by intermediate Z bosons and photons as well, 
which may also give rise to charged leptons and provide 
a background to the lepton signal from W boson decay 
when the detectors are not hermetic. Reference [17] addi¬ 
tionally derives and implements the resummation of large 
logarithms in the transverse momentum of the interme¬ 
diate W boson. 

In the calculations [17-19] the NLO corrections were 
obtained numerically in the context of a Monte-Carlo in¬ 
tegration routine. The resulting computer codes are very 
flexible in the sense that kinematic cuts on lepton or re¬ 
coil jet variables can be easily implemented. Those from 
Refs. [17] and [18] are known as RHICBOS and CHE, 
respectively, and have found wide use in comparisons to 
RHIC data. On the other hand, the Monte-Carlo inte¬ 
gration based codes are rather demanding in terms of 
CPU time. This becomes a significant drawback when 
one wants to perform a global analysis of the helicity 
distributions from the RHIC data [10, 11, 20]. Such an 
analysis typically requires many thousands of computa¬ 
tions of the spin asymmetry. Clearly, a fast and stable 
evaluation at NLO is highly desirable in this context. 

In this paper, we derive analytic expressions for the 
NLO spin-dependent partonic cross sections for elec- 
troweak boson production, including their leptonic de¬ 
cay. More precisely, we consider the cross sections di¬ 
rectly as single-inclusive lepton ones, pp —>■ where 

transverse momentum and rapidity of the charged lep¬ 
ton are observed, precisely as is the case at RHIC. We 
note that a corresponding calculation in the unpolar¬ 
ized case has been presented a long time ago [21]. We 
present a new program that produces NLO results for 
the single-spin asymmetries relevant at RHIC and out¬ 
runs the Monte-Carlo based codes by about two orders of 
magnitude in CPU time. We also include the background 
reactions involving Z bosons and photons. We expect 
our program to become a useful tool for global analy¬ 
ses of RHIC data based on Mellin-moment [10, 11, 22] 


or neural-network [20] techniques. We also use our new 
code to present comparisons of the present RHIC data to 
NLO predictions for a variety of sets of helicity parton 
distributions. 

In Sec. H we discuss the technical details of our NLO 
calculation. Section HI presents our phenomenological 
results, where we also perform comparisons with the CHE 
code of [18]. Finally, we conclude in Section IV. 

II. NEXT-TO-LEADING ORDER 
CALCULATION 

A. Eramework and outline of the NLO calculation 

We consider the single-inclusive process pp ^ ^ + X, 
where (. denotes the charged lepton (or antilepton) re¬ 
sulting from production and decay of a W boson. As dis¬ 
cussed in the Introduction, charged leptons can of course 
also be produced by an intermediate photon or Z bo¬ 
son which, subject to the experimental selection criteria, 
gives rise to a background. We hence perform all our 
calculations also for 7 and Z production and ^Z interfer¬ 
ence. For the sake of simplicity we will, however, present 
details of our calculation and explicit results only for the 
most interesting W boson case, and just highlight a few 
features specifically relevant for intermediate 7 and Z. 

We denote the momenta of the incoming protons and 
the produced charged lepton by PatPbiPi^ respectively. 
Using factorization [23], we write the polarized hadronic 
cross section dAcr which appears in the numerator of 
Eq. (1) in terms of convolution integrals of polarized and 
unpolarized parton distributions Xfa{xa, Pf), fbixb, pp) 
and the perturbative hard-scattering partonic cross sec¬ 
tions dAaab- 

dAa = ^ dXadXbAfa{Xa,PF)fb{Xb,P-F) 

a,b 

xdAaab{XaPA,XbPB,Pe,PR,lJ'F), (5) 

where 

dAaab = ^ [da~^~^ + da~^ — {da -I- da )] . (6) 

The superscripts on the right refer to parton helicities, 
so that the helicities of the second parton b are summed 
over, while he helicity difference is taken for parton a. 
The sum in Eq. (5) runs over quarks, antiquarks and the 
gluon, and the parton distributions are evaluated at the 
factorization scale pp- The partonic cross sections also 
depend on a renormalization scale ppi. The fractions of 
the parent hadrons’ momenta carried by the scattering 
partons are denoted by Xa and Xb- An analogous ex¬ 
pression for the unpolarized cross section da appearing 
in the denominator of Eq. (1) is obtained by using only 
unpolarized parton distributions and the corresponding 
unpolarized partonic cross sections, defined by averaging 
over the helicities of both incoming partons. 
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Due to the pure V — A structure of the Wqq' ver¬ 
tex, and because of conservation of quark helicity at the 
vertex, the spin-dependent partonic cross section for an 
incoming polarized quark is just the negative of the cor¬ 
responding unpolarized cross section, while for an incom¬ 
ing polarized anti-quark it is the same as the unpolarized 
one: 


dAo-qb = -ddqb {b = q',g), 

dA&qb = ddqb {h = q,g). (7) 

Note that no such relation occurs for incoming polarized 
gluons. In case of 7 and/or Z exchange, relations (7) do 
not hold. 

We now introduce the variables 


The lepton’s transverse momentum pT and its center-of- 
mass system rapidity g are related to these variables by 

y = l-^e-^ vw=^e^. (10) 

We furthermore introduce the partonic variables corre¬ 
sponding to Eqs. ( 8 ),(9): 


S=iPa+Pb), t={pa-pi)\ U={pb-piY 


t —u 

z; = 1 -F -, w= —— 
s s 1 


( 11 ) 


so that from pa = XaPA, Pb = XbPs we have 


S={Pa+PbY, T={PA-ptY, U={PB-Pif, (8) 


and 


T 


w = 


-u 

S + T ■ 


( 9 ) 


VW 1 - V 

— 5 — ~Z ■ 

VW I — V 


( 12 ) 


Writing out Eq. (5) explicitly to 0{as) in the strong cou¬ 
pling constant, we now obtain 


d^Aa 

dpxdg 


- e / f 

PT a b JvWf 

d^^ab\s,v) 

dv 


dw XaAfa{Xa, Pf) Xbfb{xb, Pf) 

i5(l — w) -\- 


a^PR) dAal2is,v,w,pF,PR) 


27T 


dvdw 


(13) 


where the dAd^J represent the leading-order (LO) con¬ 
tributions and the dAd^^J the NLO ones. 

The only LO partonic process is qq' -A- W ^ (.v^ anni¬ 
hilation, whose Eeynman diagram is shown in Fig. 1 a). 
For the NLO correction we have to include the 2 —>■ 3 real- 
gluon emission diagrams as well as the virtual corrections 
to the Born cross section. In addition, quark-gluon scat¬ 
tering contributes here as well as a new channel. Some of 
the relevant NLO Feynman diagrams are shown in Fig. 1 
(b)-(d). 

For our calculations, we work with a general (axial) 
vector structure for the IFgg'-vertex of the form 

9'^ (^9 “ «9 95) , (14) 

where Uqq' is the appropriate CKM matrix element and 
gw the fundamental weak charge. Likewise, we use a 
corresponding expression for the ITfz/^-vertex, with vec¬ 
tor and axial coefficients Vf and (and, of course, with 
Uqq' = 1). Using such general vertices will help us to keep 
better track of the couplings in the NLO calculation and 
to obtain an understanding of the underlying structure. 
Also, it allows us to extend our calculation to the case 
of 7 or Z boson exchange (for 7 Z interference one needs 


to introduce an even more general vertex structure that 
allows different couplings in the amplitude and its com¬ 
plex conjugate). The case of a W boson is recovered by 
setting Vq = Uq = \ and vi = a(, = \. 

As is very well known, various types of singularities ap¬ 
pear at intermediate stages of the NLO calculation. To 
treat these, we choose dimensional regularization with 
d = A — 2£ dimensions. This means that we have to 



FIG. 1: Feynman diagrams for heavy gauge boson production: 
a) leading-order, b) NLO virtual correction, c) NLO real emis¬ 
sion, d) NLO qg scattering. Crossed diagrams are not shown. 
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deal with subtleties that occur in Dirac traces involv¬ 
ing 75 or in the presence of the Levi-Civita tensor 
when d ^ A. 75 appears in the bDgg'-vertex (see (14)) 
and also acts as projection operator onto definite helic- 
ity states of incoming quarks or antiquarks. Likewise, 
the Levi-Civita tensor projects onto gluon helicity states. 
We adopt the ’t Hooft, Veltman, Breitenlohner, Maison 
(HVBM) scheme of [24, 25], which basically recognizes 
the four-dimensional nature of 75 and separating 

the usual four space-time dimensions from the additional 
d — A = —2e spatial ones. Technically, we compute Dirac 
traces using the Tracer package of [26]. We also fol¬ 
low Refs. [27, 28] to use a symmetrized version of the 
IT-fermion vertex. 

Because of the distinction between four- and {d — 4)- 
dimensional subspaces in the HVBM scheme, the squared 
matrix elements for the partonic processes will contain 
regular d-dimensional scalar products of the external mo¬ 
menta, but also additionally (d — 4)-dimensional ones. 
The latter have to be properly taken into account when 
the phase space integration is performed. As it turns 
out, for the unpolarized cross sections all such additional 
terms are either absent or integrate to zero, i.e. are of 
0(e) after phase space integration. However, in the po¬ 
larized case, they do contribute, and in fact a finite ad¬ 
ditional subtraction is required in the procedure of fac¬ 
torization of collinear singularities in order to maintain 
relations such as (7) beyond LO. The deeper reason for 
this is that the 75 and definitions of [24, 25], al¬ 

though algebraically consistent, cause violation of helic¬ 
ity conservation at fermion-boson vertices, which has to 
be corrected for. Since this is very well established in 
the literature (see, for example, Refs. [29-31]) we shall 
not go into any further detail here but only mention the 
salient features when they become relevant in the course 
of the calculation. 


B. Born-level cross section 


Thanks to (7), we can easily develop the calculations 
of the unpolarized and polarized cross sections in paral¬ 
lel. Up to the subtleties just mentioned, it is sufficient 
to present details only for the unpolarized case. The 
lowest-order contribution to the cross section comes from 
the 2 —)■ 2 scattering process qq' —>■ The diagram is 

shown in Fig. 1(a). As before, we use for the observed 
charged lepton, regardless of its charge. We shall see be¬ 
low that it is possible to formulate a partonic cross section 
in this generic way, despite the fact that the “lepton” can 
be either a particle or an antiparticle. We also always re¬ 
fer to the corresponding neutrino or antineutrino as the 
“neutrino” and denote it by vg. Since it remains unob¬ 
served, we integrate over its phase space. This leads to 
an overall factor 5(1 — w) for the Born cross section, so 
that 


'2^(0) 


d^d 


dd 


(0) 


dvdw 


dv 


-5(1 — w) 


(15) 


as we have anticipated in (13). Using the general vertex 
structure given in Eq. (14), we find that two combinations 
of the couplings appear in the expression for the cross 
section, which are given by 

Cl = {Vq+al){a'j + Vi) + AaqaiVqVe , 

C 2 = (Vq + aq){aj + v'j) - A aqatVqVi . (16) 

We recall that in case of an exchanged boson, we 
have Vq = Qq = V£ = = 1 and hence always Ci = 8 
and C 2 = 0. However, it is useful to keep C 2 in the 
calculation as it allows us to easily switch between W~ 
and W~^ production. The reason for this becomes clear 
when we write down the unpolarized Born cross section: 

jUqq'l^s C,(1 -v)^ + C2V^ 

dv SttN, \ V2 J (s-M^)2+r2^M^ ’ 

(17) 

where = 3, Gf = V^g^/iSM^) is the Fermi con¬ 
stant, and M\y and Tw are the W boson mass and de¬ 
cay width. Let us consider now the partonic channel 
mJ —>■ e+Fg. For this indeed Eq. (17) provides the correct 
cross section when Gi = 8 and G 2 = 0. In this way the 
cross section is proportional to ( 1 —z;)^, as required by the 
V — A structure of the interaction and angular momen¬ 
tum conservation. For du —>■ e“Pe, on the other hand, 
the cross section has to be proportional to , rather than 
(1 — u)^. This is immediately realized by interchanging 
Cl and G 2 in Eq. (17), and subsequently setting again 
Gi = 8 and G 2 = 0. Equivalently, and even more simply, 
we can just choose in (17) Gi = 8 , G 2 = 0 for ud —>■ e+Fg 
and Gi = 0 , G 2 = 8 for du —?> e~Ve to obtain the correct 
cross sections. We note that the cross sections for the 
reactions du —t e^Ve and ud —t e“Fe can be obtained by 
simple “crossing” t o zt, or u o 1 — z;. Again this may 
also be achieved by Gi O G 2 . All these considerations 
also hold at NLO, where the cross section still depends 
only on the two combinations Gi and G 2 . 

The denominator in Eq. (17) represents the standard 
Breit-Wigner form of the propagator. One often also uses 
the form (see [32]) 


(s-M^)2+r^s2/M^ ’ ' ' 

which may be obtained from the one given in (17) by the 
simple rescalings ^ M^/(H-F^/M^), Vy/Mw 
Vy/Mw! (1 + F^/M^) and multiplication of the cross 
section by 1/(1 -I- F^/M^). This also holds at NLO. 
The numerical difference between these two forms of the 
propagator is very small and negligible for our purposes. 


C. Real 2 —> 3 corrections 

At NLO, we hrst consider the 2 —>■ 3 real-gluon emis¬ 
sion process qq' —>■ where the gluon and the neu¬ 

trino remain unobserved. One of the two relevant Feyn¬ 
man diagrams is shown in Fig. 1(c). All external particles 
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can be considered as massless, so that the kinematics and 
the phase space are as usual for single-inclusive calcula¬ 
tions. The three-particle phase space in 4—2e dimensions 
may be written as [29, 31] 


dvdw 


47r 


2e 


( 47 r) 4 r(l - 2 e) V s 


X / dOi I d02 sin^ Q\ sin 62 


dz 


B{l/2,-e)Jo 




where v and w have been defined in Eq. (11) and where 
9i and O 2 are the polar and azimuthal angles of the 
neutrino in the rest frame of the neutrino-gluon pair. 
The integration variable 0 is specific to the treatment 
of 75 and in the HVBM scheme. It is given by 

z = 4fc^/(s23 sin^ 6*1 sin^ 6 * 2 ), where S 23 = su(l — w) and 
is the square of the d — 4-dimensional parts of the 
neutrino and gluon momenta, which are the same in the 
adopted frame. It is thus the only d — 4-dimensional 
invariant in the calculation [29, 31]. Note that the z- 
integral cancels against the Beta function in the last line 
of (19) for all terms in the squared matrix element that 
have no dependence on k^. 

Since the lepton pair is produced via an intermediate 
W boson, a propagator with the momentum pi + Pv^ of 
the W boson appears in the amplitude for the process. 
As a result, the squared matrix element contains 

the overall factor 


1 

(si2 - M^)2 + ’ 

with the leptons’ pair mass squared: 

S12 = {Pt+Puif' ■ 

S12 is a function of the angles 61 and 02 - Since the neu¬ 
trino is not observed, the propagator will be subject to 
integration over the phase space. We write it in the fol¬ 
lowing way: 

1 

(S12 - M^)2 + 

= ^ ^ - I - ^ 

2*5 \si 2 - - ig S 12 - -h ig 

where g = After this partial fractioning, there 

are only terms in jAI]^ with at most one power of 
S12 in the denominator, either l/(si2 — — ig) or 

l/(si2 — + ig). They are usually accompanied by 

other Mandelstam variables that also depend on di and 
02 - The ensuing terms may be readily integrated using 



( 20 ) 

( 21 ) 


the integrals 

= / d0i d 6»2 sini-^'^ 01 sin-^”^ 02 
Jo Jo 

X_ - _ (23) 

(o -I- 0 cos 01 )^ (A -I- B cos 01 -I- C sin 01 cos 02 )” 

tabulated in the Appendix of Ref. [33] . The results con¬ 
tain logarithms of various complex arguments which may 
be combined to produce manifestly real results. This 
procedure is rather tedious; we have performed numer¬ 
ous numerical checks to ensure its correctness. For terms 
with dependence on k^ the z integration in (19) is still 
trivial. The result may then be further integrated us¬ 
ing (23). The additional power of sin^ 0i sin^ 62 arising 
from the z-integral can be easily accommodated by shift¬ 
ing e —>■ e — 1 in (23). 

After integration over phase space the result for the 
real-gluon emission contribution contains singularities in 
1/e. These occur whenever we have a term in jAl]^ with 
at least a factor of 1/^3 or I/M 3 , where 

h = (Pq - Pg)"^, U 3 = iPq' - Pg)^ ■ (24) 


The poles arise when the gluon becomes collinear with 
the incoming particles, and/or when it becomes soft. The 
collinear singularities arise directly in the angular inte¬ 
grations. A soft singularity is equivalent to the invariant 
mass squared of the two unobserved particles becoming 
small, i.e. S 23 = sv(l — w) —)■ 0 , or equivalently w —)■ 1 . 
To make also the soft divergences manifest, we use the 
standard expansion 


(l-u;)-^-^ 


1 

e 


i5(l — w) -b 


1 

(!-«;)+ 


/ log(l - w) 
\ 1 — w 


+ 0{e^) , (25) 


where the “plus” distributions are defined as usual by 

[ dwf{w)[g{w)]+= [ dw[f{w) - f{l)]g{w) . (26) 

Jo Jo 

The final expression contains quadratic (1/e^) poles as 
well as single (1/e) ones. We note that due to the finite 
width Tw of the W boson, final-state singularities never 
occur. 

The NLO contributions associated with qg —)■ iiyq' 
scattering at NLO (Fig. 1(d)) can be integrated in the 
same way as described above. They develop only single 
poles in 1 /s since soft singularities are absent here. 


D. Virtual correction and factorization of collinear 
singularities 

At NLO, the interference of the virtual diagrams (see 
for example Fig. 1(b)) with the Born diagram con¬ 
tributes. As may be inferred from [27, 34] , the first-order 
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FIG. 2: Upper row: Representative initial-state collinear con¬ 
tribution for the qq' channel. Lower row: Factorization of the 
final-state collinear singularity which is necessary for the process 
with an intermediate photon. 


virtual corrections only modify the basic qq'W vertex by 
a multiplicative factor of the form 1 + 0{as). Therefore, 
when computing the interference with the Born diagram, 
the result will be twice the Born cross section multiplied 
by this factor. In our notation, we have from [34]: 


da 


(l),virt 

qq' 

dvdw 




(0),e 


dv 


-(5(1 — w) 



3 

e 


r(l-Fe) 

s ) r(l - 2 e) 



(27) 


where Cp = 4/3. It is important to take into account 
here that the Born cross section is to be computed in 
4 — 2e dimensions, where it is given by 


dd 


(0),e 


dv 


\Uqq'\^s ( GpM^ y UnV 
SnNc \ J \ s J r(l - e) 


Ciil-v)^ + C 2 V^ -C 3 S 


(28) 


Compared to the four-dimensional expression (17) a new 
combination of the vector and axial vertex factors ap¬ 
pears: 

C 3 = [ae - vif{aq - VqY . (29) 


As it turns out, this combination appears also in the 
real-emission contribution and in the factorization sub¬ 
traction discussed below, in such a way that the final 
result for the NLO correction only contains the combina¬ 
tions Cl and C 2 given in (16). We furthermore note that 
the spin-dependent Born cross section in 4 — 2£ dimen¬ 
sions with an incoming polarized quark, dAd^^q,’^/dv, is 
the negative of dAd^^^,'^/dv in (28), but with C 3 = 0 . 
This violation at order 0{e) of the relations in (7) and 
hence of helicity conservation is typical of intermediate 
results in the HVBM scheme [29, 31]. 

Adding the real and virtual contributions, the double 
poles in e cancel. We are left with single poles associated 


with collinear gluon emission. According to the factor¬ 
ization theorem, these may be absorbed into the parton 
distribution functions by a suitable subtraction which we 
perform in the MS scheme. This introduces dependence 
on a factorization scale In the upper row of Fig. 2, 
one of the two initial-state collinear situations for the 
2 —>■ 3 qq' channel is shown. Here, the variable x denotes 
the momentum fraction of the incoming quark after ra¬ 
diating a gluon. The required subtraction is of the form 
~ i Pqq®'^^^qq''‘^ 1 where Pqq is a LO Altarelli-Parisi split¬ 
ting function [35] and again the Born cross section 

for the process qq' —>■ (.vi computed in 4 — 2e dimensions. 
More precisely, in case of the contribution shown in the 
upper part of Fig. 2, in the unpolarized case, we have to 
subtract the term 

dd^°^,’^{xs,xt,u,e) 
vs dvdw Jq dv 

X Hqq{x, pp ) S{x{s 4-1) 4-u) , (30) 

where the MS scheme is defined by 


Hqq{x,pp)= - 1 -7 b- log 4? 


2 \ 


Pqqix) 


with 7 b the Euler constant and with 
Pqq{x) = Cp 


71- 

1 - 2 :+ 2 


(31) 

(32) 


Standard MS factorization requires the splitting function 
to be computed in four dimensions. After the collinear 
subtractions have been performed, we end up with the 
final NLO result in the MS scheme. 

If the incoming quark is polarized, the subtraction is 
similar, but with two crucial differences: First, one needs 
the spin-dependent Born cross section in 4 — 2e dimen¬ 
sions, given as discussed above by the negative of the 
unpolarized one in (28) but with C 3 = 0. In addition, 
as discussed in Refs. [29-31], in order to correct for vio¬ 
lation of helicity conservation in the HVBM scheme, one 
needs to use the splitting function 


APqq{x) 


Cp 


1 -F a;^ 
_(l-a;)+ 


+ 2^(1 


x) -F 4£(1 — x) 

(33) 


in the factorization subtraction. With these differences 
taken into account, the final spin-dependent NLO par- 
tonic cross sections respect the relations in (7), as they 
should. 

As already mentioned, in the case of an exchanged W 
or Z boson one does not encounter any final-state sin¬ 
gularities. Effectively, the widths of the bosons act as 
regulators here. On the other hand, for an intermedi¬ 
ate photon - which presents one of the backgrounds - 
a final-state singularity would occur if the leptons were 
massless, when the photon goes on its mass shell. Keep¬ 
ing a finite lepton mass is well beyond the scope of this 
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work and is also not necessary since the pure-photon con¬ 
tribution is in any case rather small. Also, because of 
parity conservation, it is only present in the unpolarized 
cross section and not in the single-spin one. The arti¬ 
ficial singularity that one encounters in this channel for 
massless leptons may be avoided for instance by impos¬ 
ing a cut on the invariant mass of the outgoing lepton 
pair [18], or it may be simply subtracted in, say, the MS 
scheme. Effectively, the latter approach, which we adopt 
here, introduces a (QED) photon-to-lepton fragmenta¬ 
tion function [36]. The diagrammatic situation for the 
final-state collinear splitting is shown in the lower row of 
Fig. 2. The subtraction to be performed is given by 


7 (1),photon fact 

1 

sv dvdw 


— dx 


^ , daf^^ {s,t/x,u/x,e) 


dv 


X Hij(x, fip) (5 ( s + 


t + u 


X 


,(34) 


where ddqq^^g denotes the Born-level cross section for 
the process gg —>■ 75 in d = 4 — 2e dimensions, and where 


where the coefficients A{v),B{v) may be found from 
Eq. (A.4) in the Appendix. The terms with “plus” dis¬ 
tributions represent the well-known threshold logarithms 
for the process that arise when the incoming partons have 
just sufficient energy to produce the observed final state, 
so that any substantial gluon radiation is kinematically 
inhibited. 

The other terms in the NLO result have a more compli¬ 
cated structure. The integration of terms containing (20) 
gives rise to three different types of denominators. We 
write them by introducing the function 

P(z) = . 

{zs - ' 

We then encounter the terms 

P^ = P{z,) (i = l,2,3), 

where 

Zi = 1, Z2 =w, Z3 = - - . 

1 — vw 


(37) 

(38) 

(39) 


H£y(x,fiF)=(^-^+7E-ln47r^ Pi^(x) 

with Pij(x) the appropriate 7 ^ £ splitting function. In¬ 
cluding the thus defined subtraction renders the full NLO 
cross section finite. We stress again that the pure-photon 
contribution is small, except at large lepton rapidities. It 
can in fact be vetoed experimentally because it is char¬ 
acterized by two charged leptons that almost coalesce. 
We also note that the 7 Z interference contribution does 
not produce any final-state singularities even for massless 
leptons. 

Finally, for qg scattering, there are no virtual correc¬ 
tions at 0{as). To obtain the finite cross section for 
these partonic channels, one therefore only needs the ap¬ 
propriate subtractions for the initial-state collinear sin¬ 
gularities. 



E. Final results 

Our final analytical NLO expressions for the processes 
qq' — £X, qg £X through W-boson exchange are 
presented in the Appendix. We briefly summarize a few 
features of the result for the qq' —>■ iX channel. First of 
all, it contains the usual distributions in (1 — ic), which 
dominate the cross section at w ^ 1. These multiply the 
Born cross section: 


Evidently, Pi essentially just corresponds to the propaga¬ 
tor in the Born cross section. The other two propagators 
are similar and reduce to Pi in the limit ru —)■ 1 . 

In addition to the new propagators arising at NLO, we 
also find several logarithms of the propagator terms. The 
logarithms that occur are 

, / (ws - M^)^ -I- \ 

((yEy,s-M^)^+Tl,M^\ 

+ J’ 

, f ((1-V + vw)s \ 

+T^Ma. ) • 


As seen in Eq. (A.l), they are accompanied by inverse 
tangent functions resulting from the imaginary parts of 
the arguments of the logarithms arising in phase space in¬ 
tegration. All these terms are multiplied by simple func¬ 
tions of V and w and by one of the three types of propa¬ 
gators given above. The result for the channel qg tX 
does not contain threshold distributions but does have 
logarithms of the type in Eq. (40); see the Appendix for 
further details. 


III. PHENOMENOLOGICAL RESULTS 


jAO) 

rfVP; „^i da\d, 

~ -Ojr 


dvdw 


dv 


/ log(l - w) 
V 1 — w 


4A(v) 

(l-w)+ 


+ B{v)S{l-w) , 


(36) 


We start with the unpolarized cross section for pp scat¬ 
tering at RHIC at y/S = 510 GeV. Figure 3 shows our 
LO (dashed) and NLO (solid) results for the cross sec¬ 
tion da/dpT for €+ and i~ production through interme¬ 
diate W bosons. We have integrated over |? 7 | < 1 in the 
charged lepton’s rapidity. We have used the NLO parton 



















FIG. 3: LO (dashed) and NLO (solid) cross sections at RHIC 
{y/S = 510 GeV) for G*' and ^“-production through boson 
exchange. 


distributions of [37] and the renormalization and factor¬ 
ization scales = fip = Pt- Our adopted values for 
the W mass and width are Mw = 80.398 GeV, Tw = 
2.141 GeV (later we will also use Mz = 91.187 GeV and 
Tz = 2.49 GeV for the Z boson). 

Clearly, the NLO corrections are significant every¬ 
where. They have moderate size below and around the 
Jacobian peak at pr ~ Mw/2 and become very large 
well above the peak. A close inspection of the results in 
Fig. 3 reveals a hint of a “shoulder” in the NLO cross 
sections just above px = M\y/2. This shoulder is a 
true feature of the NLO results. It comes about in two 
ways: First, the qq' channel itself has non-trivial struc¬ 
ture here. Near px = My//2, there is a complicated in¬ 
terplay between positive contributions by terms with dis¬ 
tributions in (1 — w) (“plus distributions” or J-function) 
in Eq. (A. 4), and contributions by subleading terms in 
(1 — w), among them the terms involving the functions J 
and K, which are negative around px ~ Myt/12 and be¬ 
come positive just below and above the Jacobian peak. 
This means that the q(^ channel is sensitive to the ex¬ 
act mix of positive and negative contributions. Secondly, 
the qg process makes a negative contribution below and 
around px = Mw/2. and then becomes positive. This 
intricate interplay of the various contributions is also the 
reason why the height of the peak is reduced at NLO as 
compared to LO. We note that for increasing energy y/S 
the shoulder becomes even more pronounced and in fact 
quickly turns into a double-peaked structure at NLO; see 
also [38]. This at first sight surprising feature is a mani¬ 
festation of the well-established fact [39] that the region 
around the Jacobian peak cannot be controlled within a 
fixed-order calculation. Among other things, it is sensi¬ 
tive to small transverse momenta qx of the intermediate 


W boson. There are large double-logarithmic corrections 
to the ( 7 T-distribution of W bosons at low qx, which need 
to be taken into account to all orders if one wants to ad¬ 
dress this region [40]. Such a resummation is incorpo¬ 
rated in the RHICBOS code [17]. These issues become 
relevant for precision determinations of the mass of the 
W boson from the lepton’s px spectrum near the Jaco¬ 
bian peak [41]. For RHIG, they are not really relevant 
since, if one is interested in determining polarized par- 
ton distributions, there is no need to focus on the region 
around the Jacobian peak. Rather, it is advisable to in¬ 
tegrate over a sizable range in px, so that the Jacobian 
peak region constitutes only a rather small part of the 
cross section, and to study the distribution of the charged 
lepton in rapidity. This is the strategy adopted by the 
RHIC experiments. We will therefore consider only lep¬ 
ton rapidity distributions in the remainder of this paper. 
We plan to present a more detailed analysis of the region 
around the Jacobian peak in future work. 

In order to check the validity of our analytical results 
and their numerical evaluation, we have performed ex¬ 
tensive comparisons to high-statistics runs of the NLO 
code GHE presented in Ref. [18], both for the unpolar¬ 
ized and for the polarized case. We have found excellent 
agreement. A representative example is given in Fig. 4, 
where we show the spin-dependent cross sections for £~^ 
production at RHIC, through VF+ boson exchange (left) 
and for the background channels, .Z-boson exchange and 
'jZ interference (right; the pure-photon channel does not 
contribute to the spin-dependent cross section). Both our 
analytical (solid lines) and the CHE results (histograms) 
are shown. We have followed [18] to use the polarized 
parton distributions of [11] (referred to as DSSV08) and 
the unpolarized ones of [42] which were also the base¬ 
line set in the DSSV08 global analysis. Furthermore, the 
figure is for y/S = 500 GeV, and the transverse momen¬ 
tum of the observed charged lepton has been integrated 
over the range of 20 < pt < 60 GeV. As in [18] we 
have chosen the renormalization and factorization scales 
as fiji = pp = pL = y/pp + M^/2 and assumed np = 4 
active quark flavors. In Fig. 4 the error bars of the re¬ 
sults shown for CHE correspond to numerical integra¬ 
tion uncertainties. The uncertainties in our new numer¬ 
ical calculation are smaller than the widths of the lines. 
Since our results are largely analytical whereas the code 
of [18] is based on a standard Monte-Carlo integration 
with numerical cancelation of singularities, our new code 
produces the results shown in about two orders of mag¬ 
nitude less time. Of course, Monte-Carlo based codes 
are more flexible in general, allowing the implementation 
of various additional kinematical cuts and observables if 
necessary. 

We now turn to the spin asymmetries Ap which are the 
quantities of primary interest in RHIC’s W physics pro¬ 
gram. Figure 5 shows our NLO results at VS = 510 GeV 
as functions of rj. The cross sections have been integrated 
over Px > 30 GeV, as appropriate for comparison to the 
PHENIX data [3, 4]. We have now used the new set 
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v v 

FIG. 4: Comparison of our analytical results with the corresponding ones from CHE [18] for the polarized cross sections Act for 
production through decay (left) and through intermediate Z or 7 . We have considered here pp collisions at -sfS = 500 GeV and 

have integrated over 20 < pt < 60 GeV. As in [18] the parton distributions have been chosen from Refs. [11, 42]. 


of polarized parton distributions of Ref. [10] (referred to 
as DSSV14). This set primarily contains updated infor¬ 
mation on the nucleon’s spin-dependent gluon distribu¬ 
tion, which is less relevant for weak boson production. 
However, it is also based on new results from inclusive 
and semi-inclusive lepton scattering [9], so that it offers 
new information on the quark and antiquark helicity dis¬ 
tributions as well. We use the unpolarized parton dis¬ 
tributions of [37]. The solid lines in the figure show 
our results for charged-lepton production via W decay 
for the scale choice p = M 14 // 2 , while the dotted and 
dot-dashed lines correspond to the choices p = pp and 
/i = My/, respectively. One can see that the scale de¬ 
pendence of the asymmetries is extremely weak, which 
is one of the reasons why W boson production at RHIC 
is an excellent and theoretically well-controlled probe of 
nucleon spin structure. In Fig. 5 we also investigate the 
impact of the “background” presented by Z and 7 ex¬ 
change. The dashed lines show the NLO results for the 
scale p. = Miy /2, now including the Z and photon contri¬ 
butions. As is known from previous studies [17, 18], the 
background channels dilute the spin asymmetries some¬ 
what, which is mostly due to the increase of the unpo¬ 
larized cross section in the denominator of the asymme¬ 
try. We note that the STAR experiment at RHIC is able 
to identify and subtract this background, using data as 
well as Monte-Carlo estimates, so that the data can be 
directly compared to calculations based on only inter¬ 
mediate W bosons. For comparisons to PHENIX data, 
the Zj^ background needs to be included. Figure 5 also 
shows the spin asymmetries for Z and 7 exchange alone, 
in this case integrated over 25 < px < 50 GeV corre¬ 
sponding to conditions in STAR [2]. 

Using our new NLO code, we finally compare in Fig. 6 
the results for various sets of spin-dependent parton dis¬ 
tributions to the published STAR [2] and PHENIX [3] 
spin asymmetry data taken at '/S = 510 GeV. The 


STAR AY data have been presented for various rj, sam¬ 
pled over the range 25 < pr < 50 GeV of lepton trans¬ 
verse momenta, and our theoretical results shown are 
adapted to these conditions. We note that for PHENIX 
the cut on transverse momentum is different, px > 30 
GeV, and the asymmetry is for electrons or positrons 
and hence includes the contributions from photons and 
Z bosons, as we just discussed. These are, however, rela¬ 
tively small effects (see Fig. 5), so we show the PHENIX 
data point along with our results and the STAR points. 
In view of the results shown in Fig. 5 the scale choice 
hardly matters; we use pn = pp = Myr 12,. The 
sets of spin-dependent parton distributions we use are 
from [10, 11] (DSSV08 and DSSV14), from [20] (NNPDF- 
poll.l), as well as the “statistical” parton distributions 
of [43, 44] and a much earlier set [45] known as the 
“GRSV valence scenario”. From the figure we draw the 
following observations: 

• all sets describe the W~^ asymmetry data rather well. 
The main reason for this is that the spin asymmetry is 
largely driven by the polarized up quark distribution 
which is relatively well constrained by DIS data and 
hence similar in all sets. 

• among the various sets, NNPDFpoH.l is the only one 
for which the STAR data were already included in the 
analysis, constraining the light sea quark helicity dis¬ 
tributions. As a result, the data are quite well de¬ 
scribed by the set, especially when one includes the 
corresponding uncertainty estimates [ 20 ] that we do 
not show here. Note, however, that information from 
semi-inclusive lepton scattering is not included in the 
NNPDFpoH.l set. 

• at ?7 < 0, the two DSSV sets show W~ asymmetries 
that are below the data. Since the DSSV14 set contains 
the latest information available from (semi-inclusive) 
DIS, this hints at the interesting possibility of a tension 
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FIG. 5: Single-spin asymmetries for negatively (left) and positively (right) charged leptons as functions of rapidity at v/S" = 
510 GeV. We have integrated over the range pr > 30 GeV. The solid lines show the results when the lepton originates exclusively from 
W bosons, at scale p = Mwl‘i‘- The dotted and dot-dashed lines correspond to the scale choices p = pr and p = Mw, respectively 
(note that the lines for the various scales are almost indistinguishable). For the dashed lines the background from exchanged Z 
bosons and photons has been included, using the scale p = Mwl‘i‘- Finally, the long-dashed lines show the spin asymmetries for 
Z bosons and photons alone, without the VF-boson contributions, this time for 25 < pr < 50 GeV. We have used the DSSV14 
polarized parton distributions [11] and the unpolarized ones of [37]. 


between the DIS and RHIC data, the latter favoring a 
larger Aft distribution (see also the discussion in [20]). 
It has to be emphasized, however, that we do not dis¬ 
play here any uncertainties for the DSSV set; as shown 
in [2, 11], the main DSSV08 uncertainty band is such 
that it just about touches the lower end of the error 
bars of the data points. In this sense, it is premature to 
draw any conclusions regarding such a tension. Clearly, 
it will be interesting to follow up on this issue in the 
context of a new global analysis, especially when addi¬ 
tional experimental information becomes available. 

• in the framework of the statistical parton distributions, 
the helicity distributions are obtained along with the 
unpolarized ones and depend on only very few param¬ 
eters to be determined from data. As one can see from 
Fig. 6 (and as discussed in [44]), the model describes 
the RHIC data quite well. 

• the GRSV valence scenario of [45] describes the W~ 
asymmetry data strikingly well. The main distinctive 
features for this set are assumptions about the breaking 
of SU(3) in the relations between nucleon spin structure 
and hyperon /3-decays, and the ansatz [46] 

Au(a;,(5g) Ad(a;,(5§) 

at a low initial scale Qq. Since Au and Ad are known 
to have opposite sign, the latter ansatz forces the ra¬ 
tio Ad/Au to be negative. This requirement, along 
with the condition Au -F Ad < 0 imposed by the DIS 
data and the assumptions about SU(3)-breaking, is re¬ 
alized in this model by a fairly large positive Au dis¬ 
tribution and a negative (and even larger in absolute 


value) Ad one. Evidently, the STAR data prefer such 
a sizable positive Au. We note that one of the sets of 
Ref. [47] has a similar Au distribution and hence de¬ 
scribes the W~ asymmetry data similarly well [18]. It 
will be interesting to see whether also the large negative 
Ad of [45] is realized; unfortunately, the Ad contribu¬ 
tion to the W~^ asymmetry is typically overwhelmed 
by the Au one. Note that a negative Ad pulls the W~^ 
asymmetry to more negative values (see (4) in the in¬ 
troduction), which may explain why the GRSV valence 
scenario shows the most negative asymmetry of all the 
sets at 77 < 0. Needless to say that the GRSV valence 
scenario has not been confronted with the latest (semi- 
inclusive) DIS data. 


IV. CONCLUSIONS 

We have presented a new analytical NLO calculation 
of the partonic cross sections for single-inclusive lepton 
production at RHIC, when the lepton originates from 
the decay of an intermediate electroweak boson, espe¬ 
cially a W boson. Our numerical code based on ana¬ 
lytical phase space integration is much faster than exist¬ 
ing Monte-Carlo integration based codes. In this way, 
we hope that our code will be a valuable tool for future 
global analyses of the proton’s helicity parton distribu¬ 
tions that include the new high-precision data for A'^ 
asymmetries obtained at RHIC. Our results may also be 
useful to obtain insights into the analytical structure of 
the partonic cross sections, for example in terms of their 
threshold logarithms or their behavior in the vicinity of 
the Jacobian peak. 
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FIG. 6: Comparisons of NLO results for AY for various sets 
of helicity parton distributions [10, 11, 20, 43, 45] to the STAR 
data [2] taken at v/S = 510 GeV and to the PHENIX mid¬ 
rapidity points for electrons/positrons with \ri\ < 0.35 [3]. The 
cut 25 < pt < 50 GeV has been applied on the lepton's trans¬ 
verse momentum. Note that the PHENIX points are for pr > 30 
GeV and includes the contributions from photons and Z bosons. 
We have chosen the scales pa = pp = Mw/‘2. 


We have also presented new comparisons of the latest 
RHIC data with the NLO predictions for some of the sets 
of polarized parton distributions available in the litera¬ 
ture. In line with observations in the earlier literature we 
have found that the data prefer a rather sizable positive 
Au helicity distribution in the proton. 
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Appendix 

In this Appendix, we present some of our explicit NLO 
results. We first consider the qq' channel when an inter¬ 
mediate W~ boson is produced (for example through du 
scattering), for which effectively Ci = 0, (72 = 8 in (16) 
(see discussion after Eq. (17)). We define the functions 


K{z) 

J{z) 


arctan f 1 ^ ~ 

\zs — J 


= log 


_i_ 


2 M\y 


w 


K(z) 


(A.l) 


with the usual (Heaviside) step function. In addition to 
the values zi = 1 ,Z 2 = w,Z 3 = (1 — v)/(l — vw) of 
Eq. (39), we introduce 


zo = 0, Zi = 1 — V + vw 


and we set 


Ji — J — Al(2i) . 
We then find for production of a W~: 


(A.2) 


(A.3) 


^ KVP f GpM^ 

dvdw t:Nc \ >/2 




2{1 + wY 


2^ ( l0g(l -W)\ ^ Pqq(w) 


V 1 — w 


- 2 log(l - vw) 

_l_ ^ F 


+ + Q + 21og(l - w)^ log <^(1 - li') + {Jo - J 2 - J 3 + Ji + k{Ko - K 2 - K 3 + Ki)) 

V f Jq — 2 J 3 + J 4 Jo ~ Ji 


2 \ 1 — vw 

1 1 A w^ 


1 — V A vw 


+ V^{P 2 


\ I-w J , Up K^S J 


2 1-w 


(^Jo ~ 2^2 A J 4 + ~ ^ 4 )^ } 1 


u —f 1 — vw, w 


1 — V 
1 — vw 


(A.4) 


with the splitting function Pqq of Eq. (32), and with Note that despite appearance the expression is perfectly 

well regularized at rc = 1. 

^ _ 2Mw{Yw + ^w)_ 
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By applying crossing one obtains the corresponding 
cross section for q'q —)• W~g. Crossing is achieved by 
changing v ^ 1 — vw, ic —)■ (1 — n)/(l — vw) and multi¬ 
plying the result by the Jacobian v/{l — vw). We do not 
give the crossed result explicitly here. 

Writing the NLO partonic qq' cross section for general 
Cl and C 2 in the form 

Cl + C 2 ddf ^ , (A.6) 

we find that crossed- Since the result for 

production is obtained in our calculations by setting 
Cl = 8 ,C 2 = 0 (see Sec. II B), we thus have 

= Sddi^^^ = dd-^] , 

qq'^W+g 1 qq'g ’ 


dd 


( 1 ) 

q'q^W+g 


dd{ 


( 1 ) 


= dd 


crossed 




We remind the reader that the cross section for a 
polarized incoming quark differs just by a sign from the 
corresponding unpolarized one (see Eq. (7)) while that for 
an incoming polarized antiquark involves no sign change. 
The cross sections for intermediate Z bosons may be con¬ 
structed from (A. 6 ), using (A.4) and its crossed variant 
and inserting the appropriate coupling factors Ci and C 2 
in each case. 

Secondly, we also present the result for the channel 
gq —?> W~(f in the unpolarized and the polarized case: 


s_d^^ 

dvdw 


Tr I Uqq' I ^ f GpM^ \ 2 


ttA^c V 


V^P2 < 2 (1 - w)w - Pqg{w) 


Ji) — 2 J 2 + Ji -\ -(ATo — 2 K 2 + K 4 ) + 2 log ( . ^ 

w V^(I ~ 


-2 


s d^Ad, 
dvdw 


gg 


TR\Uqq'\^ / GpM^ ^ ^2 

ttA^c V V 2 


P 2 2 (I - w) - APgg{w) 


Jo — 2 J 2 + d 4 -I {Kq — 2 K 2 + K 4 ) -\- 2 log ( —— 

W \^(1 ~ w)s 


r 


(A. 8 ) 


where Tr = 1/2 and 

Pqg{x) 


^ + (1 - , 


We note that the terms in square brackets have a similar 
structure as the penultimate one in (A.4). Finally, for 
qg W~q' we find 


APqg{x) — 


2 ( 22 ^-!) • 


(A.9) 


sdJd^qg 

dvdw 


TR\Uqq>\^ ( CfM^ 

V ^2 


2 r 


V r 2M^ 

^ (Jo - 2 J3 + J4 + k{Ko - 2K3 + Ki)) 


1 — vw 


2„,,2 


-I- P3V W 


2(1 — w)w — Pqg (ui) ( Jq — 2 J 3 -I-J 4 -I—^ {Kq — 2 K 3 774 ) + 2 log 


fJ'F 


;(1 — w)i 


I — V — vw + 2v^w — v^w^ (1 -|- nw)(l — 2n -I- vw) 


M 2 


— (Jo — 2 J 3 -|- J 4 ) 

- 1 - V 

(1 — -I- nw)2 

-I- -(Jo — d 4 )(l — w)(l — 2v + 2vw) — w(l — V — 2w + vw) 
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1 — vw 


(Jo — J4 + r(Ko ~ K4)) 


1 — vw 

1 — 3u -I- 2^2 -I- — Sv'^W + v'^w'^ 
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(A.IO) 
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where 


and 


The corresponding spin-dependent cross section for an 
incoming polarized quark again just differs by a sign; 
(A.ll) see (7). 


K = 


r2 Sf2 
-1^ VI/ i- 


(A.12) 
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